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, ' We suggest a new mean-field method for studying the thermodynamic competition between mag- 

^ _ netic and superconducting phases in a two-dimensional square lattice. A partition function is con- 

■ structed by writing microscopic interactions that describe the exchange of density and spin fluctua- 
' tions. A block structure dictated by spin, time-reversal, and bipartite symmetries is imposed on the 

_^ ' single-particle Hamiltonian. The detailed dynamics of the interactions are neglected and replaced 

by a normal distribution of random matrix elements. The resulting partition function can be calcu- 

1 i ■ lated exactly. The thermodynamic potential has a structure which depends only on the spectrum of 

quasiparticles propagating in fixed condensation fields, with coupling constants that can be related 
directly to the variances of the microscopic processes. The resulting phase diagram reveals a fixed 
^ ■ number of phase topologies whose realizations depend on a single coupling parameter ratio, a. Most 

5-H ' phase topologies are realized for a broad range of values of a and can thus be considered robust 

with respect to moderate variations in the detailed description of the underlying interactions. 

2 . 

^ ■ PACS numbers: 71.10.Fd,71.27.+a,74.25.Dw 

S ■ I. INTRODUCTION 

Studies of high-temperature superconductors have revealed a rich phase diagram, with coexisting magnetic and 
superconducting correlations. These phase structures can be complex and include d-wavc pairing, stripes, or the 
pscudogap phenomenon. [l|, 0, H, E| Theoretical models and numerical studies on a lattice indicate that the richness 
of the phase structure results from a delicate energy balance between competing states. @, [|| It follows that model 
predictions can be sensitive to small changes in the parameters of the theory or to details of the numerical approach. 
The question then arises of which properties of the phase diagram are constrained by the basic underlying symmetries 
and which are sensitive to the detailed dynamics of the interactions and to numerical approximations. 

The purpose of this paper is to address this question with a new mean-field approach. The method is based on 
random matrix theory and consists in constructing a Hamiltonian that retains the basic spin, time-reversal, and 
bipartite symmetries of the problem but simplifies the dynamics of the interactions considerably. The theory is 
radically different from the familiar Hubbard or t-J model. Here, we construct the model at a deeper microscopic 
qq ' level and describe interactions that are mediated by the exchange of density and spin fluctuations. This construction 
is inspired by random matrix models of the strong interaction, for which the QCD interactions are mediated by 
I* ■ single-gluon exchange. Although natural in QCD, a microscopic description involving bosonic fields may be more 
. 5^ | controversial in the context of high-T c superconductors. Such an approach is similar to low-energy effective theories of 
^ i antiferromagnets and superconductors (til. M. M. l9l. [Tpj or to the antiferromagnetic spin-fluctuation exchange theory, pdl 
EH EH, 0] Here, in contrast to these models, no particular assumption is made regarding the detailed dynamics of 
the exchange fields. Instead, we adopt a coarse description in which the block structure of the interaction matrix is 
dictated by the underlying symmetries of the Hamiltonian while individual matrix elements are drawn at random. 

A random matrix approach offers three advantages. First, since the theory is constructed at a more microscopic 
level, it allows us to relate global properties of the phase diagram to specific microscopic mechanisms. Second, the 
simplified dynamics produces a mean-field model that can be solved exactly: the gap equations are polynomial. Their 
roots can therefore be studied as a function of the coupling parameters of the theory. Third, in the vicinity of 
critical points, the thermodynamic potential has a Landau-Ginzburg form in which the expansion coefficients satisfy 
specific symmetry constraints inherited from the deeper microscopic level. These constraints help us to identify those 
topologies that can be realized in the system and rule out those that violate the constraints. 

The motivation for applying methods used in QCD to the high-T c problem follows from the strong analogies existing 
between these systems. First, the restoration of chiral symmetry with increasing quark density can be understood from 
the analogous behaviors of QCD and metamagnets. The chiral condensate plays the role of a staggered magnetization 
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which vanishes abruptly as an external magnetic field is increased, driving the system through a first-order phase 
transition. fl5j Second, single-gluon exchange is attractive in the antitriplet channel and leads to the Cooper pairing 
of quarks. This form of pairing can lead to a long-range order called color superconductivity. [Tol. IT7L fisj] The degrees 
of freedom that are involved in pairing are different from those involved in the chiral broken phase of QCD, so that 
color superconductivity competes directly with the breaking of chiral symmetry. 

Random matrix methods have been extensively applied for studying the phase diagram of QCD. [TBI. fl9l. l20l. [2TI. [2^| 
In recent works . [23l l24l , I25I [26j| we studied the phase diagram of QCD with three colors and two flavors as a function of 
temperature and quark chemical potential. The partition function was constructed as an integral over random matrices 
that mimicked the basic structure of quark-quark interactions but neglected their detailed dynamics. These matrices 
were given a block structure that reflected the spin, color, and flavor symmetries of one-gluon exchange. Inside a given 
block matrix, no further correlations were assumed among the matrix elements, which were then drawn at random 
on a normal distribution. This approach produced a mean-field model that could be solved exactly. The resulting 
effective potential gave polynomial gap equations whose roots could be determined analytically or numerically. The 
effective potential contained a single free parameter, defined as a coupling-constant ratio that measured the relative 
strength of the interactions in the chiral and diquark channels. As this coupling ratio was varied, the phase structure 
passed through a restricted number of distinct topologies. Moreover, starting with coupling constants with the values 
appropriate for single-gluon exchange, large variations were required to alter the topology of the phase diagram. We 
thus concluded that the QCD phase diagram was robust against moderate variations in the detailed dynamics of the 
interactions. 

In general, random matrix models are useful in providing a global picture of the phase diagram. Being mean- 
field in nature, such a picture is only a starting point that requires the proper inclusion of thermal, quantum, and 
spatial fluctuations (which play an important role in the high-T c problem as a consequence of the Mcrmin- Wagner 
theorem[27j) if it is to be quantitatively reliable. Nevertheless, the random matrix approach can be useful in providing 
an overview of the strength of the order parameters and their sensitivity to coupling parameters and can thus serve 
as a means for identifying those characteristics of the phase diagram that are protected by symmetry. 

It is worth noting that the "random" part of the theory has nothing to do with disorder. Instead, the philosophy here 
consists in constructing a Hamiltonian with a block structure dictated by the underlying symmetries of the problem 
and replacing individual matrix elements by random variables, conventionally drawn on a normal distribution. This 
construction can be regarded as equivalent to integrating over many Hamiltonians that meet fundamental symmetry 
requirements but which differ from one another in the detailed implementation of the dynamics of the interaction. 
The theory can thus appear elaborate at first sight since it involves a large number of statistical variances. However, 
the number of free parameters decreases at each step of the calculations so that the theory becomes simpler as one 
proceeds towards the solution of the problem. In fact, the final form of the thermodynamic potential depends only 
on a single parameter ratio, and its functional form has a well-defined structure that can be interpreted in terms of 
quasiparticle energies. Hence, most of the effort involved in constructing the model is "upfront" but worthwhile since 
it provides relationships between the global phase diagram and the microscopic parameters of the theory. Overall, 
the procedure is relatively simple and could be implemented in many other problems. 

We will consider a fermion system on a two-dimensional square lattice and construct its partition function at finite 
temperature and finite chemical potential. In Sec. [TTJ we show that extension of the methods used in QCD poses some 
challenges. First, while the basic interactions of QCD arc naturally formulated in terms of quarks exchanging gluons, 
no such natural description is available at an elementary level for the cuprates. We will thus explicitly assume that 
the interactions can be described as the exchange of density and spin fluctuations. The structure of the interaction 
between the fermions and the fluctuation fields is dictated by S , ?7(2)-spin, time-reversal, and bipartite symmetries. 
We will show that this formulation leads to a four-fermion effective potential whose terms can be compared to those 
of the Hubbard model. A second challenge is the need to account for the d-wave character of the superconducting 
order parameter, which forces us to introduce an explicit momentum dependence of the fermion states. As we are 
seeking to construct a model that suppresses as many details of the interaction as possible, we limit ourselves to a 
"coarse-grained" momentum description in which the first Brillouin zone is divided into four sectors in order to mimic 
the symmetry patterns of the antiferromagnetic and superconducting order parameters. 

Section IIIII is devoted to the derivation of the effective potential. We find that it is possible to construct a theory 
in which antiferromagnetism and d-wave superconductivity are favored and compete, while the s-wave channel is 
repulsive. Such a theory gives greater statistical weight to spin-fluctuation fields with a large momentum exchange. 
When deriving the effective potential, we will see that the theory simplifies considerably as one proceeds through 
the calculations. The initial model of Sec. |TT] involves as many as eight different variance parameters; the final 
thermodynamic potential depends on a single parameter ratio. Its interpretation in terms of quasiparticle energies is 
also given in Sec. IIIII 

We study the resulting phase diagrams in Sec. [iVl where we show that, as was the case for QCD, there are only a 
finite number of possible topologies. As in the QCD problem, the topology of the phase structure changes gradually 
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as the ratio of coupling constants is varied. Section [V] contains a summary of our main findings and our conclusions. 

II. CONSTRUCTION OF THE RANDOM MATRIX MODEL 

We consider a system of electrons on a two-dimensional square lattice and model the competition between magnetic 
and superconducting orders as a function of the chemical potential, ^t, and the temperature, T. We construct 
interactions that satisfy spin, time-reversal, and bipartite symmetries with simplified but integrable dynamics. This 
approach ensures that the properties of the model arise solely as a consequence of the symmetries. 

A. Order parameters and the parameterization of momenta and frequencies 

We wish to define random matrix correlators that mimic the basic structure of antiferromagnetism and supercon- 
ducting order parameters. Working at finite temperature T in an imaginary time formalism, the antiferromagnetic 
order parameter assumes the form 

m AF=( V'i(P + Q,W Jl )or ai 3^(P ! W n ) V (1) 

\pw„a/3 / 

where p are momenta in the first Brillouin zone, Q = (nh/a, irh/a) is the AF ordering vector (a is the lattice spacing), 
cu n = (2n + 1)ttT are fermion Matsubara frequencies, a and (3 are spin indices, <x are the spin Pauli matrices, and 
(...)= Tr(. . . e~^ H ) denotes a thermal average. Similarly, the d-wave order parameter is given as 



msc-d 



9(p)^(P,^n)ipihP, -w„) ) , (2) 
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where 



is the d-wave form factor. 

Our aim is to construct correlators that mimic the momentum couplings in Eqs. (TIJ and j2| ona coarse level. First, 
we divide the Brillouin zone into four regions related to one another by either a momentum shift in Q, 

{p}^Q{p} = {p + Q}- (4) 

or by momentum reversal, 

{p}^P{p} = -{p}. (5) 

Next, we replace the exact form factor g in Eq. ([2|) by the simplified form factor 

<f> d (p) = sign(#(p)) , (6) 

where 4>d is a crude approximation of g. It neglects its variation with momentum, which is related to the detailed shape 
of the wave function, but exhibits the same d-wave symmetry as g, changing sign for every ninety degree rotation in 
the Brillouin zone. We believe that such an approximation captures the essential symmetry of the problem and is 
sufficient to describe (i-wave pairing. 

The particular form for <p d guides our parametrization of momentum states. Four momentum regions are chosen 
as the sectors in which <p d has a given sign. One possible division of the Brillouin zone is shown in Fig. [TJ The 
approximate form factor <pd is —1 in regions 1 and 3 and +1 in regions 2 and 4. States in regions 1 and 3 are related 
to those in regions 2 and 4 by a shift of ±Q. Regions 1 and 2 are related to regions 3 and 4 by momentum reversal. 
Inside each region, states are labeled by an index i = 1, . . . , M, where M scales with the total number of lattice 
sites and M — * oo in the thermodynamical limit. These states also count the different Matsubara frequencies. From 
one region to another, the states are parametrized as follows. If a given index i refers to a state (p, LJ n ) in region 
1, then the states labeled with the index i in 2, 3, and 4, respectively, correspond to (p + Q, u> n ), (—p,—ui n ), and 
(— p — Q, — u> n ). (Note the change of sign in the frequencies of the last two terms.) 
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FIG. 1: Parametrization of the first Brillouin zone. Regions 1 and 2 are related by a momentum shift by ±Q = ±(nh/a, nh/a), 
as are regions 3 and 4. Regions 1 and 3 and regions 2 and 4 are related by momentum reversal. 



With this parametrization, the AF order parameter is written as 

m AF = ( H ^l^cLfi {Paf) t ,„ Sn 1> Bdtfl J , (7) 

\r,a=lt,j=Xa,j9=T,| / 

where r and s are region indices and the four- by-four matrix Taf couples momenta separated by Q: 

10 

T AF = (a 1 ) Q ®(l) P = I J o o i I • (8) 

10. 

Similarly, the superconducting order parameter has the form 

/ 4 M V 

m S c~d = ( X] zJ ( T sc-d) r , s Sij ip 3tj> i J (9) 

\r,s — 1 — 1 / 

where Tsc-d is now blind to shifts by Q and couples states with opposite momenta with a sign dictated by the 
d— wave form factor of Eq. ©: 

0-10 

(°-i)p = I _°i o o I 
10 0. 

Note that with our parametrization, the Kronecker symbol Sij in Eqs. (|7|) and © automatically selects equal Mat- 
subara frequencies for the AF order parameter and opposite frequencies for the superconducting SC order parameter. 

Similar arguments can be followed for an s— wave order parameter. With an isotropic form factor </> s (p) = 1, the 
corresponding momentum projector is given as 



Tsc-s = (1)q®(cti) 



p 



( 1 
1 
10 

Vo i o o 



(11) 



and the s-wave order parameter has the form given by Eq. (0) with Tsc-d replaced by T$c- 



5 



B. Constraints imposed by symmetries 



We now turn to the construction of the random matrix interactions. Inspired by random matrix models for 
QCD[23|, [2H, we write the partition function of the system as a path integral 

Z(fi,T) = J V^Vi> dH int P(H int ) e -^(Ho+H lM )i> ) (12) 

where ip' and tp are independent fcrmion fields, i?; n t is a matrix describing the random interaction with elements 
distributed according to the distribution P{H m t) to be defined below, and Ha is the non- interacting part of the 
single-particle Hamiltonian which contains temperature and chemical-potential terms. 

This formulation is radically different from the more familiar Hubbard or t-J models, where the interaction terms 
are described by effective four-fermion potentials. Here, instead, the interactions are described at a deeper microscopic 
level so that fermions interact with fluctuation fields via current terms, ip^H^tp. This formulation is directly inspired 
by the Yang-Mills Lagrangian for the strong interaction mediated by gluon exchange. In the context of the high-T c 
problem, these fluctuation fields can describe either interactions carried by phonons, antiferromagnetic fluctuations, 
or more complex effective interactions. Here, we will not attempt to identify the nature of these fields nor to specify 
their dynamics. We will assume only that interactions can be described in terms of such fields and consider how their 
description is constrained by the symmetries of the system. 

Interaction terms. According to the parametrization introduced in Sec. Ill Al the fermion fields are described 
by spinors with eight components (four momentum regions and two spins). Similarly, the random matrices iJi n t are 
composed of 8 x 8 block matrices of size M x M. Due to the symmetries of the system, these blocks are not completely 
independent. 

Consider first the constraints imposed by SU(2). For a fixed pair of momentum region indices, e.g., (r, s), we can 
write 



-f^int, r,s,f,T -^int, r,s,f , J. 
Hint, r,s, -Hint,r,s, J., J. 



^a^H^, (13) 



where 07, = (l,er), H Kr s with /i = represents density-fluctuation fields, and H KT}S with /1 = 1,2,3 describes 
spin-fluctuation fields. The interaction Hamiltonian can be made invariant under a spin unitary transformation, U, 
by requiring that the vector (Hi- rtS , H2, r ,s, -ff3;r,s) simultaneously undergoes a space rotation i?, with R satisfying 
U^OiU = Rij<jj. The partition function itself, Eq. (|12[) . is then made invariant under SU(2) transformations by 
requiring that the probability distribution, P{H mt ), is invariant under the corresponding spatial rotations R. 

Consider next constraints related to time-reversal invariance. Following the work of Monthoux on spin-fluctuation 
exchanges [l3j . the fields are taken to be real in coordinate representation — they do not carry an electric charge. 
As will be seen below, the integration over these fields produces a four-fermion interaction that contains squares of 
density terms, ~ {tp^ip) 2 ! an d spin currents, ~ {tp^ ertp) 2 , thus leading to a two-body potential with a time-reversal 
symmetry. Now, because the fluctuation fields are real, the Fourier components of their matrix elements must 
satisfy the constraint 

(p f w n \Hp\q,u m ) = (-p, -w„|i/ M | -q, -u m )* ■ (14) 

If one divides the Brillouin zone into two subspaces of states {p,Lu n } (regions 1 and 2) and {— p, — lo„} (regions 3 and 
4) and adopt the parametrization introduced in Sec. Ill A[ this condition can be cast in the form 

= ph;p, (15) 

where 

p-(?S) («> 

reverses both momentum and frequency. Hence, the matrices must have the block-structure 

Hfj, I ^* r>* j , (17) 

where are Hermitian and C M are complex symmetric. (40| 
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Finally, we turn to the bipartite symmetry, appropriate for a square lattice, composed of two interpenetrating 
sublattices A and B. Consider the transformation 

, , , / +V>(r) if r G A, . . 

#0-(_^ r ; ifreR (18) 

The full Hamiltonian H is not expected to be invariant under such a transformation because the kinetic terms couple 
fields defined on neighboring sites. However, we will assume that the interaction part of the Hamiltonian, H ult , 
is bipartite invariant (as is the case for the U term in the Hubbard model). We first determine the momentum 
representation of Eq. (| 18|) . Using the definition of Q = (71-, 7r) 7i/a, Eq. (fT8|) can be rewritten as tp(r) 1— > exp(iQ • 
v/h) ip{r). Then, in the momentum representation, the bipartite transformation takes the form 

^(p)^Q^(p) = V(p + Q). (19) 



Dividing the Brillouin zone in the two subspaces {p} (regions 1 and 3) and {p + Q} (regions 2 and 4), the bipartite 

= QH^Q, (20) 



invariance of is written as 



with 



Hfj, must then have the block-structure 



1 j ; i 



where and are Hcrmitian. An alternative, but weaker, requirement on could be that matrix elements 
between states with momenta pi + Q and P2 + Q arc equal to those between pi and p2- This amounts to take 
equal diagonal blocks in the right side of Eq. (|22|) . with no constraint on the off-diagonal blocks. This second choice 
results in a free energy with a slightly different form but leads to the same main results as the choice of Eq. 
This alternative form is discussed further in the Appendix. 

Combining the requirements imposed by the three symmetries, we arrive at an interaction matrix of the form 



Hint = CT fi. 

fi=0 



^ B^d B^ lo C^d C^ \ 

Bfio B^d C^ C^d 

S~1* O* D* 

U (io D nd D fio 

, s~l* sy* r>* o* 

\ /!0 £>fj,d J 



(23) 



where the 4x4 block structure of iJi nt refers to the regions of the Brillouin zone. Each block is described by an 
M x M matrix; blocks B^d and i? MO are Hcrmitian, while C^d and C MO are complex symmetric. Overall, H mi contains 
16 independent blocks. 

Probability distribution. We represent the 16 blocks by Ab with b = 1, . . . , 16. Their matrix elements are drawn 
on a normal distribution 

P(H int ) = exp (-8M S ' Tr ( A b4) ) . ( 24 ) 

where represent inverse variances. This form allows us to perform the integration over H lnt analytically and thus to 
determine the partition function exactly. [4l| In order to make the partition function invariant under SU(2) rotations, 
the inverse variances associated with each of the three blocks that describe a spin fluctuation exchange are chosen 
equal. Since 12 of the 16 independent blocks describe spin-fluctuations and four describe density fluctuations, we 
arrive at a total of 4 + 4 = 8 independent variances. The resulting distribution function is then given as 



P(H int ) = exp(-8Af(E| 0<I '&(S M B^) + S| d '&(B (i .Bi) 

+E| 0o Tr(i? 0oJ Bt o ) + S| o Tr(B ■ Bj) + E^Tr^C^) + £^Tr(C d • C+) 

+ Z 2 COo Tr(C 0o Cl) + S^ o Tr(C • Cj))) , (25) 
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where the inverse variances can be tuned individually at will in order to favor various scattering mechanisms. 
Non-interacting terms. The non-interacting part of the single particle Hamiltonian is written as 

^H Q ^jj = V> f (-/i + (It + T t ) ip • (26) 

Here, the chemical-potential term, /i, is a scalar while (It and T t are matrices which describe temperature dependence 
and hopping terms, respectively. 

Temperature is introduced via Matsubara frequencies. Following previous work in QCD.[23l. l2~i| we include only 
the two lowest frequencies, ±«7rT. With this approximation, (It takes the form 

T = diag(i7rT, -iirT) ® (l) spin ® (1) Q ® (a 3 )p, (27) 

where the first term on the right side is an M x M diagonal matrix. Hence, in each momentum region, half the 
states have a positive frequency while the other half have a negative frequency. The final term on the right side of 
Eq. (|27|) serves to implement the parametrization introduced above: states with a fixed label i in regions 1 and 2 have 
frequencies opposite to the corresponding states in regions 3 and 4. 

Although limiting the sum over Matsubara frequencies leads to an oversimplified description of temperature depen- 
dence, we believe it to be sufficient to determine the general characteristics of the phase transition. In fact, we will see 
below that the parameter T serves as an energy scale which influences the energy balance between the various order 
parameters, much as the average thermal energy does when all frequencies are taken into account. The inclusion of all 
frequencies would certainly modify the resulting phase diagram but only in the trivial sense that every temperature 
T is mapped monotonically to a new value. [28j We will not seek to refine the T-dependence here since such mapping 
does not alter the phase topology, i.e., the occurrence and the order of transition lines. Since T as introduced here is 
an arbitrary temperature scale, we will also drop the factor tt to simplify notation. 

The hopping term T t in Eq. Q26[) is written as 

r t = diag(t, -t) ® (l) S pin ® (ojOq ® (l)p, (28) 

where the first term on the right side is a diagonal M x M matrix containing M/2 elements of +t and M/2 elements 
of — t. The matrix T t mimics the nearest-neighbor hopping energy 

£ p = -2to (cos(p x a/?i) + cos(p y a/h)) , (29) 

whose band is symmetric around £o = and which satisfies £ p +q = — £ p . Here, again, we neglect the detailed 
momentum dependence of the hopping energy and retain only its symmetries around £o = and under a shift by Q. 
In Eq. (|28|) . t is a measure of the strength of the hopping term, which can be tuned against the variances and thus 
against the interaction strength. 

III. THERMODYNAMIC POTENTIAL 

In this section we evaluate the thermodynamic potential corresponding to the partition function of Eq. (|12p using 
methods which are standard in random matrix theory. [3, [23| It is worth emphasing that the random matrix model of 
Eq. (12) is solvable, i.e., the saddle-point method to be discussed below becomes exact in the thermodynamic limit. 
We will concentrate on the main results and present more detailed calculations in the Appendix. These calculations are 
performed in three steps. The first step consists of integrating over the interaction matrix elements according to the 
distribution P(H int ) of Eq. (|25p. This integration leads to a four-fermion interaction, Y, which receives contributions 
from the eight independent blocks of H m t in Eq. (TJ3|). In the second step, the fcrmion fields arc rearranged and 
the terms of Y are put in two groups. The first group contains products of bilinears of the form ~ ViV', (with an 
implicit sum over i), which includes the antifcrromagnctic order parameter. The second group contains products 
involving bilinears of the form i/>iipi, which are relevant for superconductivity. To keep track of momentum and spin 
indices while rearranging terms, we use the generalized Fierz identities of the Appendix. In the final step, quartic 
fermion terms are linearized by means of Hubbard-Stratonovitch transformations, which introduce an auxiliary real 
field a, to be associated with an antifcrromagnctic order parameter, and a complex field A, to be associated with 
superconductivity. The resulting form of the partition function is given as 



Z{ft,T)= / dadAdA*e- 8Mnia > A \ 



(30) 
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where n is the thermodynamic potential 

fi(<r,A) = A|A| 2 + Bo 2 - - log((^ 2 +t 2 - m) 2 + |A| 2 + T 2 ) 



-^^((v^+^ + M^ + IAp + T 2 ), (31) 



with 



A - 8 ^2 ^2 I" ^2 v2 ^ !" v2 ^2 h f^ - ' ( 32 ) 

X^BOd ^BOo ^COd ^COo L B ^C d ) 

( i i i i i i i i V 

S = 8 ~^2 ^2 ^2 ^2— + y2~ + + + ' ( 33 ) 

\ ^BOd ^BOo ^COd ^COo ^Bj Zj B a ^Cj L C„ / 

Equations (f3"Tj) - ([3"3|) constitute the main result of the model. Note that the factor 8M in the argument of the 
exponential in Eq. (|3T)|) plays the role of the volume of the system since M scales with the total number of lattice 
sites. The phase diagram can be established with the help of the saddle-point approximation to Eq. (piO)) . which 
requires the simultaneous solution of the two gap equations dfl/da = and dfl/dA = 0. Due to the logarithm term 
in the right side of Eq. (|31[) . these equations are polynomial in the auxiliary fields and can be solved analytically. 
In the thermodynamic limit M — > oo, the saddle-point approximation becomes exact and the solutions of the gap 
equations that achieve the lowest value for f2 describe the thermodynamic phases of the system. 

As mentioned in the introduction, the model becomes increasingly simple as one proceeds through the calculation. 
From the initially large number of parameters required to describe H m t in Eq. (|23|) t there remain only a few constants 
in the final form of the thermodynamic potential [Eq. (|31[) ]. This potential has a remarkably simple structure. 
Anticipating the results of the following section, we note first that the topology of the phase diagram depends only 
on a single parameter ratio 

a = -, (34) 

whose strength characterizes the relative importance of superconductivity and antiferromagnetism. The terms ^4|A| 2 
and Ba 2 represent the energy cost of creating a constant field in the corresponding channel. As shown in the Appendix, 
the logarithmic term in f2 corresponds to the determinant of the Hamiltonian for a single fermion in fixed constant 
external fields a and A. These terms have the generic form 

~^log(iT- £± )HT- £± ), (35) 



where e± are the quasiparticle energies 



2 



1/2 



e±= (Vt 2 +a 2 ± t i) +|A| 2 • (36) 



These expressions are strongly reminiscent of the quasiparticle energies of earlier mean- field models. |29l. l30| Neglecting 
the triplet order and approximating the square of the form factor as 4> 2 d 1 in Eqs. (10) and (11) of Ref.[29|], one finds 

quasiparticle energies of the form E±(p) = ((y^ + a 2 ± fi) 2 + | Ap) 1 / 2 with a = 2Jm and A = Jd. Equation (|36|) 

shows a result of similar structure with, however, the simplification £ p ~ t which is a consequence of our coarse 
description of momentum states which leads us to ignore the detailed momentum dependence of kinetic-energy terms. 

Thus, the basic structure of the potential f2 is simply related to the energies of the elementary excitations of the 
system for fixed constant external fields. This means that the thermodynamic potential could have been constructed 
immediately from the knowledge of the quasiparticle energies alone, without going through the steps described in 
Sec. In] Note, however, that the additional information which comes from constructing the interactions at the more 
microscopic level is useful. Through the dependence of A and B on the individual variances, it establishes connections 
between the microscopic mechanisms and the global properties of the system. 



IV. PHASE DIAGRAM 



We now consider the various topologies that can be realized in the phase diagram. Despite the simplicity of the 
thermodynamic potential, a full exploration of the parameter space is a considerable task. Thus, we concentrate on 
a restricted number of physically relevant cases. 
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First, we assume that the interactions are attractive in the antiferromagnctic channel. As shown in the Appendix, 
this requires that the variances satisfy the inequality 

B = 8 (-^-^-^-^ + ^ + ^ + ^ + yT- I >°> ( 37 ) 

\ ^BOd ^BOo ^COd ^COo ^B d L B„ ^C d / 

which implies that spin-fluctuation exchange is stronger than density-fluctuation exchange. In fact, this condition 
can be related to the requirement that the on-site potential in the Hubbard model is repulsive. In momentum 
representation, the Hubbard Hamiltonian is written as 

H K = ^2zpiptpip ap + U <Wp3,P2+P4 V , | Pi V , Tp2^Ip3^1p4' ( 38 ) 

"p P1P2P3P4 

with U positive. The corresponding partition function contains the thermal average Tr(e~ HH/T . . .), where the equiv- 
alent of the weighting factor e~ Hii / T in the random matrix model is the term e Y . Hence, a positive U in the Hubbard 
model would correspond to a negative four-fermion term of the form 1 l ; ]i' l l'y' l p\j'4'ii m the random matrix model. Close 
inspection of the four-fermion potentials in Eqs. (|A3|) - (|Al0j) reveals that they contain such terms and that they will 
be negative provided the inequality of Eq. (|57|) is satisfied. 

As a second assumption, we consider interactions which are attractive in the superconducting d-wave channel but 
repulsive in the s-wave channel. This requires a particular choice of the inverse variances. In the Appendix, we show 
that such an interaction can be found among those that favor spin-fluctuations by putting more statistical weight on 
the blocks B G and C D than on B,j and Cd- This choice leads to a positive constant A in Eq. (|3"!2)) . Such a choice of 
the variances is expected to result from interactions which favor large momentum transfer of order ~ Q. This result 
can be related to the antiferromagnetic spin-fluctuation model of Ref.fbij in which the spin susceptibility is peaked 
at a momentum exchange ~ Q. In the present approach, however, we do not attempt to describe the dynamics of 
spin exchange in detail but rather use the coarse device of the inverse variances to tune the relative strengths of the 
exchange mechanisms. 

Given the restricted parameter space that results from these two assumptions, the system can develop different 
phase structures in the (/i,T) plane as a function of the parameter ratio a = B/A. The phase structures can 
be grouped according to their topology. We identify four distinct topologies which emerge as the parameter a is 
gradually increased. The system switches from one topology to the next at specific values of a that depend on the 
strength of the hopping term t. In general, systems with larger t develop a larger kinetic energy per charge carrier 
and are found to favor superconductivity over antiferromagnetism. 

It is useful to note that because our model only includes nearest-neighbor hopping, it cannot distinguish between 
hole- and electron-doped systems. This symmetry can be seen in the potential of Eq. (|31[) . where f2 is an even function 
of fi. Hence, /! = corresponds here to half-filling. In principle, next-nearest neighbor and higher-order hopping terms 
could be added to the model by using a hopping matrix, T t , whose eigenvalues have a sign distribution that reproduces 
the symmetries of the corresponding kinetic energies in the various regions of the Brillouin zone. Although we have 
limited ourselves to nearest-neighbor terms here, we expect that a model with a more elaborate hopping matrix would 
distinguish between electron and hole dopings. 



A. Antiferromagnetism alone 

For the smallest values of a = B/A, superconductivity is too weak to compete with antiferromagnetism. Such a 
situation corresponds, for instance, to a large value for A, leading to a prohibitively large energy cost ~ A|A| 2 for 
creating a constant field A. 

Consider first the system at zero temperature and at half-filling (fi = 0). Setting A = in Eq. (|3"Tj) . we find the 
gap equation 



= 2Ba - = 0, (39) 

A=0,/j=0,T=0 a > 1 



which gives either a = (paramagnetic phase, PA), or a = \Jl/(2B) — t 2 (antiferromagnetic phase, AF). The latter 
solution is real if the hopping term is not too strong, 2Bt 2 < I. Given that B scales as an inverse variance [see 
Eq. (|33p]. this condition is equivalent to the inequality t/txH < 1> where the threshold value txH = I/V2B is a 
measure of the interaction strength. When this condition is fulfilled, a = y / l/(2B) — t 2 is the absolute minimum of 
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CI and the ground state is antifcrromagnetic. In the rest of this work, we will explicitly assume that t < ixH, so that 
the half-filled state is antiferromagnetic. 

In contrast to our approach, mean-field studies of the Hubbard model find an antiferromagnetic ground state no 
matter how weak the interactions are or how large the hopping term is.[3ll.l32T| There, the absence of a threshold results 
from the logarithmic singularity that occurs in the density of states at the edge of the magnetic Brillouin zone. At half- 
filling, the divergent density of states leads to a gap with an exponential dependence on t/U, a ~ t exp(— I'KyJTJTj^), 
where U is the strength of the on-site repulsion. In the random matrix approach, the detailed band structure is 
ignored and, because of the coarse description of momentum states, no divergence appears in the density of states. As 
a result, the interaction must be sufficiently strong to produce an antiferromagnetic ground state. Given the level of 
approximations underlying our approach, this behavior is not unreasonable. Discrepancies with microscopic theories 
are to be expected in cases where the condensates are weak and thus sensitive to fluctuations. Understanding the fate 
of these condensates clearly requires more than mean-field approximation. 




FIG. 2: Phase diagram for t/tTH = 0.5 and a < 0.1. The transition from the antiferromagnetic phase (AF) to the paramagnetic 
phase (PA) is second-order at half-filling (thin line) and first-order at zero temperature (thick line). These two lines merge at 
a tricritical point, t. Here, temperature is plotted in units of T c = T c (fj, = 0, t = 0), which is the transition temperature at 
half-filling in the limit t — ► 0, while chemical potential is plotted in units of ao = l/v2B, which represents the AF field at 
half-filling, zero-temperature, and for t — * 0. 



We now turn to the phase diagram in the (fi, T) plane. Choosing i/i-p/H = 0.5, the superconducting phase does not 
develop so long as a < a a with a: a «0.1. The corresponding phase structure is shown in Fig. [5] and resembles that 
of chiral symmetry breaking in QCD with two flavors and three colors, in the limit where a color-superconducting 
phase is ignored (see Ref . |23l] ) . The gap equation dQ/da = 0| A=0 has a form similar to the QCD problem. It has a 
solution with a = 0, describing a paramagnetic phase, and solutions which satisfy the quadratic equation 

x 2 + 2x(-^ +T 2 -^) + (^ + T 2 ) 2 + = 0, (40) 

where x = a 2 + t 2 . For moderate fj, and high temperature, the only real solution is a — and the system is in a 
paramagnetic phase. Decreasing T at fixed n, two additional real solutions 

a = ±a AF = ±L 2 -t 2 -T 2 + ^ + ^Vl-64BVT 2 ) ^ (41) 



can be found below the critical temperature 



1 



T c (fx, t)=[ — -^-t 2 + —^\- 16iV + 64B V <2 ) • (42) 



1/2 
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Below T c (fx, t), these solutions are local minima of the free energy fi(cr, A = 0) and a = becomes a maximum. The 
finite roots ±o~af describe an antiferromagnetic phase. They vanish at T c (/x, i), which thus characterizes a second- 
order transition from an antiferromagnetic to a paramagnetic phase. In the opposite regime of low temperatures and 
finiteu, the transition is found to be discontinuous. It takes place along a first-order line [actually, a triple line (see 
Ref.|23()] which starts on the zero-temperature axis at 



/il = \TiT + i2 ' (43) 

and extends with decreasing \i toward the zero-/x axis. This line meets the second-order line T = T c (p,,t) at the 
tricricital point t = (/«3,T3), given as 

1/2 



.3 = -^f^^^l , (44) 



8B 8B 



1/2 



The first- and second-order lines meet with equal slopes dT/dfx. 

The thermodynamic potential 0(cr, A = 0) can be expanded as a series of powers in a near the critical lines. 
The result resembles a Ginzburg-Landau expansion. For < fj,^ and near T c (fx,t), the free energy is found to 
scale as f2(<7, A = 0) « f2(cr = 0, A = 0) + bi(fi, T)a A + 0(a e ), so that the second-order phase transition has the 
critical exponents of a mean-field (f> 4 theory. Near the tricritical point, one has fi(c, A = 0) ~ fl(<r = 0, A = 
0) + bs([i,T) a 6 + 0(a 8 ) and the critical exponents are now those of a mean-field <fi 6 theory. Note, however, that the 
coefficients 64 and be here are known functions of /i and T. 



B. Competition between antiferromagnetism and superconductivity 



T/T c 
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FIG. 3: Phase diagram for t/t-rn = 0.5 and a = 0.2. In addition to the phase structure of Fig. [2] there is a superconducting 
phase (SC) emerging out of the antiferromagnetic phase via a first-order transition. This phase undergoes a second-order 
transition to the paramagnetic phase at either higher fi or higher T. Here, t is a tricritical point. The scales T c and 00 are 
those defined in the caption of Fig. [5] 



As the coupling ratio B/A increases, a superconducting phase can be favored over an antiferromagnetic one. As 
doping is increased, the new carriers disrupt the antiferromagnetic correlations. If the interaction is sufficiently strong 
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in the pairing channel, this can lead to a transition to a superconducting phase with A / 0. Such a transition can 
generally take place in one of two ways, either through the appearance of a "wedge" of mixed broken symmetry 
with both iT^O and A ^ 0, with continuous transitions toward the pure AF and SC phases, or via a discontinuous 
transition between the two pure phases. The first case was encountered in a random matrix model of Ref. [23| when 
the coupling ratio of QCD was altered in favor of color superconductivity. Here, however, the second case is found as 
shown in Fig. [3l 

To understand the onset of superconductivity, consider a pure phase with a = and A ^ 0. The gap equation, 
dVl{<j = 0, A)/g?A = always has a root A = 0, which is a local minimum at either large /i or large T. For moderate 
/i and for T less than 



the solution A = becomes a local maximum, while O(0, A) exhibits two local minima, given by the roots with 



which describe the superconducting phase. A$c vanishes on the curve T = T c /\(n,t), which is thus a second-order 
transition line. In particular, the thermodynamic potential behaves like Q(a = 0, A) « 51(0, 0) + cla\ A| 4 + 0{\ A| 6 ) in 
the vicinity of the phase boundary. 

The curve T = T cA (fx,t) meets the T = axis at ^+ = (1/(4,4) + t 2 + a/I + 16At 2 / \AA)) 1 / 2 . For t 2 < 1/(2A), 
this line also meets the fi — axis, so that the region T < T c A.(n,t) contains the half-filled state with \i = T = 0. 
For t 2 > 1/(2A), on the other hand, the curve T = T cA (/i,t) has the shape of a dome, which starts on the T = 
axis at /z_ = (t 2 + 1/(4A) - a/I + IQAt 2 / (AA)) 1 / 2 , reaches a maximum T max = (1/(8A*))a/I + l6At 2 at ^i max = 
a/64A 2 < 4 — 1/ (BAt) and decreases again to meet the T = axis at \x = 

The condition for the superconducting phase to emerge out of the AF phase is now clear. The curve T = T c a must 
end at fi + > fj,i, where /ii is given by Eq. (|43p . In that case, there is an intermediate region where the superconducting 
state achieves a lower energy than both the paramagnetic and the antiferromagnetic states. For a fixed ratio i/frrHi 
the threshold condition /i + = /ii gives the critical value a = oja which marks the onset of superconductivity. The 
result is a decreasing function of t. with cka ~ 0.28 for t = 0, cka ~ 0.1 for t = 0.5 txH, an d o;a ~ 0.03 in the limit 
t — > ixH- For values of a > oa, the superconducting phase develops in a wedge adjacent to the antiferromagnetic 
phase, with a first-order transition toward the antiferromagnetic state at the lower [i and a second-order transition 
toward the paramagnetic phase at the higher /j,, as illustrated in Fig. [3J 

Figure [4] shows the zero-temperature auxiliary fields oaf and A$c as a function of the chemical potential, /z. 
Note that A$c vanishes at fj, = )A+ with the mean-field exponent 1/2 so that A$c ~ (/•*+ ~ m) 1 ^ 2 - The variation 
of the antiferromagnetic field ct^f with /i should be not be considered significant; it is a direct consequence of the 
approximate description of temperature dependence. In fact, taking the sum over all Matsubara frequencies produces 
a constant condensation field by a mechanism similar to that observed in the phase diagram of QCD with two colors 



As a = B/A increases above a a, the superconducting phase boundary slides up along the first-order transition line 
between the antiferromagnetic and the paramagnetic phases. At a new critical value, a = ctb, the superconducting 
phase boundary reaches the tricritical point. The value of is readily determined from the condition T c a(^3, t) = T3 
where fx^, T3, and T c a are, respectively, given by Eqs. (|4"4"1) . (|4"S"|) . and P5)l . For t — > 0, this gives ctb — v / 2/2; at, then 
decreases with t, is equal to a& ~ 0.62 for t = 0.5 fTH and reaches ctb ~ 0.29 in the limit t — > £th- 

For a > at, the topology of the phase structure is changed. The two second-order lines T c (n,t), Eq. (|4"2")l . and 
T c a(a*, i), Eq. (|46p . now intersect at a new critical point, 6, with 




(46) 




(47) 



and light masses (see Ref . [28j ) . 



C. Bicritical point 




(48) 



(49) 
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FIG. 4: Zero-temperature auxiliary fields as a function of fi for the parameters corresponding to the phase diagram of Fig. [3] 
The scale <tq is defined in the caption of Fig. [2] 



T/T c 




FIG. 5: Phase diagram for t/tTH = 0.5 and a — 0.8. The superconducting phase has grown beyond the tricritical point of 
Fig. [3] The two second-order lines (thin lines) separating, respectively, the antiferromagnetic and superconducting phases from 
the paramagnetic phase now meet at a bicritical point, b, which is also the end point of a first-order line (thick line) between 
the AF and SC phases. The scales T c and ao are those defined in the caption of Fig. [2] 



In the vicinity of b, the thermodynamic potential can be expanded as f2(cr, A) r* 51(0, 0) + a 4 |A| 4 + 64 cr 4 + C4 er 2 |A| 2 , 
where 04, 64, and C4 are known coefficients satisfying the inequality 40464 — c| < 0. As a result, the global minimum of 
£1 can only be realized by a pure phase cither paramagnetic, antiferromagnetic, or superconducting. A mixed-broken 
symmetry state with both non-vanishing a and A cannot be a global minimum. Therefore, (/ih,Th) is a bicritical 
point which ends a first-order line separating the antiferromagnetic and superconducting phases. The coupled gap 
equations actually have a root with a mixed-broken symmetry. We have verified that, away from b and at lower 
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temperatures, (i) this state is always metastable (i.e., it is a saddle point) whenever its fields are real and (ii) the 
transition line between the antiferromagnctic and the superconducting phases is first-order all the way to the T = 
axis. The resulting phase structure is shown in Fig. [5] for the case t = 0.5 £th and a = 0.8. 

For even larger values of a, the bicritical point migrates toward the /i = axis and reaches it when a = 1. In 
this case, Tt, = T c (p = 0, £). The antiferromagnctic phase thus exists only at half-filling, and the superconducting 
phase develops for all finite /i in the region T < T c &([i,t). For a > 1, the superconducting phase wins over the 
antiferromagnctic state, which disappears from the phase diagram. 

D. Discussion 

Symmetries of Hi nt and correlations. The interaction part of the Hamiltonian was constructed by imposing 
three different symmetries, which provide the correlations necessary for generating finite condensates. Spin rotational 
and time-reversal symmetries are intimately related to magnetism. (39l| The spin symmetry is obviously necessary in 
order to express the order parameters. Time- reversal symmetry leads to a block structure which produces terms of 
the form iplipjiplipj, thereby yielding pairing condensates ~ iftjipj and ipl^l- The bipartite block-structure induces 
correlations among states whose momenta are separated by Q. The resulting four-fermion potential, ~ ipjipj^ipi, 

thus yields condensates of the form ~ where the coupled momenta are separated by Q. 

Effect of the hopping term. Since the form of the thermodynamic potential is relatively simple, it is easy to 
vary the hopping parameter in order to understand its influence on the phase structure of the system. Figures [6] — [8] 
show the phase diagrams that are realized by a system with t/tTH = 0.7 and a = 0.05, 0.2, and 0.8, respectively. 
The resulting phase diagrams exhibit the same basic topologies found in the previous case. No additional phase 
diagrams are introduced by a larger value of t. Increasing the hopping term has two primary effects. First, the 
threshold parameters «a and on, are both reduced. Second, larger t tends to increase the region occupied by the 
superconducting phase and decrease that occupied by the antiferromagnetic phase. This can be understood directly 
from the form of the quasiparticle energies s± in Eq. (|36|) . In the antiferromagnctic phase, e± directly depends 
on the combination a 1 + t 2 . Increasing t without excessively increasing e± thus requires a reduction in a. In the 
superconducting phase, however, e± varies with t as (t ± /i) 2 + | A| 2 so that the /x term tends to reduce the sensitivity 
of A with respect to variations of t. 

0.8 i 1 1 1 1 1 1 1 1 
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FIG. 6: Phase diagram for t/frrH = 0.7 and a = 0.05. The transition from the antiferromagnetic phase (AF) to the paramagnetic 
phase (PA) is second-order at half-filling (thin line) and first-order at zero temperature. These two lines merge at a tricritical 
point, t. The scales T c and <ro are those defined in Fig. [2] 

Absence of a mixed-broken symmetry phase. In Sec. IIV C[ we found that fl has a series expansion of the 
form fl(a, A) w O(0, 0) + a4|A| 4 + 641T 4 + C4(J 2 |A| 2 near the bicritical point, where the coefficients 04, 64, and C4 
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FIG. 7: Phase diagram for t/tTU = 0.7 and a = 0.2. Thin lines are second-order and the thick line is first-order; t is a tricritical 
point. The scales T c and ctq are those defined in Fig. [2] 




FIG. 8: Phase diagram for t/t-rH = 0.7 and a = 0.8. Thin lines are second-order and the thick line is first-order. & is a bicritical 
point. The scales T c and <ro are those defined in Fig. [2] 



are known functions of t and a which satisfy the inequality 40464 — c| < 0. As a result, a phase with both non- 
zero a and non-zero A is never realized. Such a mixed-broken symmetry state is in fact a solution of the coupled 
gap equations dil/da — and dfl/dA = but is unstable since the matrix of second derivatives has a negative 
determinant (i.e., the mixed-broken symmetry phase is necessarily a saddle point of £l(a,A)). In this regard, the 
present random matrix result differs from mean- field results for both the Hubbard [2^, [33| an d t~J models, [3(| for 
which mixed-broken symmetry phases are found along with a tetracritical — rather than a bicritical — point in the 
T) plane. Coexis ting phases have also been reported in a number of numerical studies including variational cluster 
perturbation theory. [34f and variational Monte Carlo methods, [3^] while other studies report phase separation. (36j 
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Calculations from cellular dynamical mean-field theory show a strong tendency to a homogeneous coexistence phase 
at weak coupling and a first-order phase transition at strong coupling. [37j 

The differences arising at mean-held level between the random matrix approach and microscopic models may well 
be a consequence of the neglect of density of states effects in the potential of Eq. (f5Tj) . In the absence of a logarithmic 
divergence in the gap equations, condensation fields are weakened and the energy balance between the phases can 
be upset. We note, however, that some authors have questioned whether the phase diagram is necessarily controlled 
solely by the Van Hove singularity in the bare density of states [Hj as mean-field microscopic results would seem to 
show. The discussion of this point shows an example of results that should not be considered universal or robust; 
further careful numerical studies are needed to settle this issue. Although a random matrix model will not give a 
definite answer, its comparison to other mean-field approaches can help in identifying those features that are sensitive 
to the specifics of models and their numerical treatment and which are thus not protected by the symmetries. 

Absence of an exact higher symmetry. It is interesting to ask whether the potential of Eq. (pTlj) can exhibit 
a higher symmetry, such as the 5*0(5) symmetry proposed by Zhang, [38[ either exactly or approximately in certain 
regions of the phase diagram. Such a symmetry should manifest itself in the possibility of writing O as a function of 
a single combination of the condensation fields, such as a 2 + |A| 2 . 

If we restrict our attention to a = B/A < 1 so that AF order truly competes with SC order, the answer is negative. 
Near the bicritical point, we found that the coefficients in the series expansion of the thermodynamical potential, 
O ph 0(0, 0) + 04! A| 4 + fo4<j 4 + C4I A| 2 <7 2 , satisfy the inequality 40464 — c| < 0. This inequality eliminates the possibility 
that the fourth-order terms are a perfect square. Therefore, no simple symmetry — allowing for a rotation of the AF 
order parameter into the SC order parameter — can be identified at the level of two-body correlators. Moreover, we 
have been unable to identify any other symmetry from an expansion of O to higher order. 

The limit a — > 1 can roughly be seen as one of approximate higher symmetry. In fact, the bicritical point migrates 
to the vertical axis, fib — > 0, and we have 4a4&4 — c| cx /i 4 — > 0. Then, can be approximated as!l si O(0, 0) + 
(-y/aX I A| 2 + \/bl(7 2 ) 2 with 04 ~ 64, a form which shows an approximate higher symmetry In the limit a — > 1, the 
antifcrromagnetic phase exists only in the immediate vicinity of half-filling and the superconducting phase dominants 
the region of finite chemical potential. Realizing such a situation in actual materials would require very specific 
relationships between the coupling constants. 

V. CONCLUSIONS 

In this paper, we have suggested a mean-field model for investigating the thermodynamic competition between 
magnetic and superconducting orders in a two-dimensional square lattice. This model describes interactions at 
a more microscopic level than the familiar Hubbard or t-J Hamiltonian through the introduction of density and 
spin-fluctuation exchanges. The single-particle Hamiltonian is given a block structure that is dictated by spin, time- 
reversal, and bipartite symmetries, and the detailed dynamics of interactions are replaced by a normal distribution 
of random matrix elements. The model is formulated in a momentum representation so that a coarse description of 
the first Brillouin region allows us to introduce a d-wave form factor that possesses the appropriate sign symmetry 
but neglects detailed momentum dependences. Although this approach may seem elaborate at first sight, the model 
simplifies dramatically as one proceeds through the derivation. In fact, the resulting thermodynamic potential has a 
simple and well-defined structure that depends solely on the form of the quasiparticle energies in given condensation 
fields. 

We have explored a number of physically relevant cases for which the interactions are attractive in both the 
antifcrromagnetic and the d-wave channels and repulsive in the s-wave channel. Such interactions naturally place 
greater weight on spin-fluctuation exchanges, particularly those involving a momentum transfer ~ Q. The phase 
diagram is found to depend on a single parameter ratio, a, and a limited number of topologies appear as a function 
of a. None of these topologies allows for a mixed-broken symmetry phase with coexisting antifcrromagnctism and 
superconductivity, probably as a consequence of the absence of the singularity in the bare density of states. Except 
for the smallest values of a, which result in a phase diagram involving antiferromagnetism alone, the range of values 
for a that correspond to a given topology is rather large. As a result, these phase topologies should be regarded as 
robust with respect to moderate variations of the detailed description of the interactions. 

The random matrix approach described here has a broad range of applications, as we have demonstrated by 
studying the phase diagrams of QCD and that of the cuprates. Given the relatively simple structure that is obtained 
for the thermodynamic potential, this method is convenient for obtaining a direct zeroth-ordcr description of the 
phase structure. In fact, such a potential could also be written immediately given only knowledge of the quasiparticle 
energies. However, the explicit construction of the microscopic interactions carries additional and possibly useful 
information regarding the dependence of the coupling parameters on microscopic variances. In the vicinity of critical 
points, the thermodynamical potential can be expanded as a power series of the condensation fields in a form similar to 
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a Ginsburg-Landau theory. The information relating microscopic processes to the coupling parameters thus imposes 
significant constraints on the Ginzburg-Landau coefficients and offers an improved understanding of the relation 
between global properties of the system and its microscopic description. 

In some sense, this paper can be regarded as "open source theoretical physics." It contains an "algorithm" that can 
be adapted for use in other problems. The shortcoming if this approach is clearly that it is no more than mean-field 
theory. Its merit lies in the fact that it gives results which are averaged over an ensemble of theories. Model-dependent 
details arc thus eliminated. The resulting thermodynamic potential is dictated by the underlying symmetries of the 
problem and is likely to be robust. 

Here, we have concentrated on the possible topologies of the phase diagram and hence on the global minimum of 
O as a function of /j, and T. We believe that the results presented here represent a generic mean-field phase diagram. 
Other features found either in theory or experiment are likely to be "fragile." For instance, the absence of a mixed- 
broken symmetry state in the present model suggests that this state is sensitive to mo del- dependent details, numerical 
approximations, or to the detailed properties of the sample studied. Indeed, the mixed-broken symmetry state is not 
generic in the cuprates.[37l| More comprehensive studies of f2 can be useful in revealing, e.g., the competition between 
local minima in f2 and other relatively fragile structures that could be sensitive to model-dependent details and 
numerical approximations. In this sense, we suggest that the present methods may provide a useful complement to 
the investigation of detailed models of these systems. 

APPENDIX A: CALCULATION OF THE THERMODYNAMIC POTENTIAL 

In this appendix, we sketch the calculations leading to the thermodynamic potential, £1, of Eq. (|31[) . The calculations 
closely follow those of the random matrix model for QCD given in Ref. [23[ . 

We start with the single-fcrmion Hamiltonian of Eq. (|2"3")) and integrate exp(— ip> H- m tip) over the matrix elements 
of H lnt according to the normal probability distribution P(Hi nt ) of Eq. (|25f . This integration can be written as 

J DH int P(H int ) e -^ H ^ = e Y , (Al) 

where the four-fermion potential, Y, receives a contribution from each of the eight independent block matrices 

Y = Y B0d + F Bd + Y BOo + F Bo + Yood + Y Cd + Y COo + Y Co . (A2) 
An integration over the matrix elements describing density fluctuations gives 

r-L— J2 (!* h + 2, 2, + 3, 3, + 4, 4. t ) (h.c.) , (A3) 



Yi 



BOd 



32MEi 



Y, 



32A/^|7 X) ^ 2 * + 2 * ^ + ^ 4 * + ^ 3 *) ( h ' C -) ' ( A4 ) 

Yc0d = 32MS 2 ? (li 3j + 2l 4j + lj 3l + 2j 40 (h - C ° ' (A5) 
^— J2 (!* 4 J + 2 * 3j + lj 4, + 2j 3,) (h.c.) , (A6) 



32MS^ 0o ( 



where (h.c.) indicates Hermitian conjugation and li lj is a compact notation for ^ a V'liaV'ija- (Here, ipij a represents 
a state with the momentum label j in region 1 and a spin index a.) Integration over the matrix elements of the blocks 
describing spin fluctuations produces the terms 

YBd = 32ME 2 S ( U a ^ + 2 ' a 2 > + 3 ^ n 3 * + ^ a ^ ' ^- C ') ' ^ A7 ) 

B d y 

32A l^ 2 E (*< a 2 i + 2* a l i + 3 J a 4 « + A J a 3 ') ' ( h - c -) ' ( A8 ) 

B„ y 

E a ^ + 2 > a ^ + h a 3 ' + 2 J a ^ ■ ( h ' c -) ' ( A9 ) 

Cd ij 

0* °" 4 j + 2 t a 3j + lj <x 4, + 2j a 3,) • (h.c.) , (AlO) 



Yb = 



Y c 



32ME n 

y 
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where the notation lj <x 2j stands for ^2 a g ipiia <T af3ip2j/3- 

Next, fermion fields are rearranged to make the condensation channels apparent. Schematically, terms of the form 
ipjipjipjipj (with the sum over i and j implied) are brought into the form — Tpj^Pj, which gives rise to condensates 
~ ipjipi that are relevant for antiferromagnetism. Similarly, terms of the form ipfyjiplipj arc rewritten as 

which contain condensates of the form ~ ip i and are relevant for superconductivity. Uncrossing necessitates keeping 
track of momentum and spin indices. To uncross spin indices, we use the SU (2) Fierz identities 

5 ab S cd = l -8 ad 5 cb + - ( r ad - ( T cb , (All) 
3 1 

CafCcd = -^S a d5 c b — - <T ad • <T cb (A12) 

for the antiferromagnetism channel and the identities 

SabScd = - S ac 6db + 2 <Tac ' adb > (A13) 
11 3 1 

CTafCcd = g fiacSdb + ^ l) ac{v l) db ~ ^ ( CT 2 )ac(c r 2 )rfh + 7J i* 7 ?) ac{& 2.) db (A14) 

for superconductivity channels. To uncross momentum indices, each quadratic term ipj,ip s (where r and s now denote 
momentum indices) is written as an element of a 4 x 4 momentum matrix which is then decomposed on a complete 
basis of 16 Hermitian Dirac matrices [23]] that satisfy 

Tr(r fc r z ) = 44; (M = 16). (A15) 

There exists a representation in which this basis contains the operators T A p, ^SC-d and Tgc-s, which were introduced, 
respectively, in Eqs. (|8|), (fT0|) . and (|11| . as well as an operator that is proportional to the hopping term, r t , of Eq. (|28|) . 
With this complete basis defined, we uncross momentum according to the Fierz identities 

(Fk)ab(Tl)cd = / ^ y Xklmn(^m)ad(^n)cb, (A16) 
mn 

Xklmn = Tr (r fc r„r;r m ) (A17) 

lb 

for terms relevant for antiferromagnetism and 

(Tk)ab(Tl)cd = %klmn(Fm)ac(Fn)dbi (A-18) 

Xkimn = T7. Tr (r fc r„rf r m ) (Ai9) 

lb 

for terms relevant for superconductivity. Here, Tr denotes a trace and Y T is the transpose of T. 

Restricting ourselves to the antiferromagnetic and superconducting channels, uncrossing gives us the four-fermion 
interaction 

Y = Y AF + Ysc-d + ■ ■ ■ (A20) 

= a A F (^T AF a^j) 2 + a S c- s {^Tsc-sV-z^) (V> Tsc-s^i 1 ) 

+asc-d (-0 t rsc-dcr 2 V ,t ) ("0 Tsc-d^ip ) , (A21) 

where a summation over spin, momentum, and random matrix indices is implied in each fermion bilinear. The 
coefficients in the respective channels are given by 

a - 1 ( 1 1 1 + 1 + 1 + (A22) 

256M y ^sod ^BOo ^COd ^COo ^B d ^B ^C / 

1 / 1 + 1 + 1 + 1 3 3 3 3 \ 

256M y^ 2 B0d S| 0o ^C0d S C0o S B d ^c d ^c / 

1 ( 1 1 1 1 3 3 3 3 \ , A „„ 

asc-d = Trr--7 ^2 ^2— + ^2 ^2 ^5- + ^5- - ^o- + ^5- • ( A24 ) 
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Here, a positive (negative) coefficient corresponds to an attractive (repulsive) channel. From the above equations, we 
thus see that interactions for which the inverse variances S 2 ^, S^, £c d , or are smau compared to ^ 2 B0d , T,% Qo , 
^cod' or ^COo are attractive in the antiferromagnetic channel. Such interactions favor the exchange of spin fluctuations 
over the exchange of density fluctuations. Similarly, interactions can be made attractive in the d-wave channel and 
repulsive in the s-wave channel by favoring the elements of the off-diagonal block matrices B„ and C Q over those of 
the diagonal blocks and CV Such interactions favor the exchange of spin fluctuations with a large momentum 
transfer, ~ Q. 

Below, we will assume that the random interactions favor antiferromagnetism and superconductivity in the d-wave 
channel and are repulsive for s-wave pairs; this situation corresponds to a particular choice for the variances such 
that aAF > and asc-d > whereas asc-s < 0. We will thus neglect the s-wave channel in the remainder of the 
calculations. In this case, combining Eqs. (fT2"]) and (|Al[) yields a partition function of the form 

Z(jx,T) = J V^VtP e'^H^+Y^ (A25) 

where, according to Eqs. (|A20|) and (IA21I) . 

Y ~ a A F {ip^TAFa^) 2 + a S c-d (ip^ S c-d<J2^) (ip ^sc-d^ip ) ■ (A26) 

The quartic fermion terms can now be written as the difference of two squares. Each square is linearized by the use 
of a Hubbard-Stratonovitch transformation, 

e AQ2 ~ J dzcxp (~- Qx^j , (A27) 

which introduces an auxiliary field x. When applied to Eq. (|A26[) . such transformations introduce a real field, a, 
to be associated with antiferromagnetism and a complex field, A, to be related to superconductivity. The partition 
function is then written as 



Z(jm,T) 



[ dadAdA* [ V^Vtp exp (-^H ip — 

J J V ^asc-d ^qaf 



X 



/ /\ A* \ 
exp I -a^TAF^ - —ip^T S c-d^2W ^sc-d^i' J , (A28) 



where Ho = —\x + f2y + T t [see Eq. (f2l))) ]. Using the spinor W = (tp^ , ipt), the fermion bilincars can be arranged in the 
Gorgov form fy^Hty, with 

H=( -^ + ^T + T t + aT A F - iAT sa \ (A29) 

y iA*T S c-d fi-U T -T t + (tTaf ) 

Then, integrating over the fermion fields gives 

V^ViPe"^^ = Det[H], (A30) 
so that the partition function in Eq. (|A28j) can be written as 

Z(fx, T) = J dadAdA* e -8A/n( ff ,A) ) ( A31 ) 
where f2 is the thermodynamic potential 

Q = ^7 (~^- + — - 4 log Bet(H )) . (A32) 
32M \a SC -d a A F J 

Calculating the determinant of H in Eq. (|A29[) then gives 

0(cr,A) = A|A| 2 + Ba 2 - i log(( v / ( r 2 + 1 2 - M ) 2 + |A| 2 + T 2 ) 

-ilog(( v ^ 2 T^ + M ) 2 + |A| 2 + r 2 ), (A33) 
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with 



^ = 8 ( — - — + ^- ^- ^ + ^- ^ + ^H . ( A34 ) 



V^i 

N^BOd ^BOo ^COd ^COo ^B d ^Bo ^C a / 

n ( 1 1 1 1 1 1 1 1 \ 1 /AN 

5 = 8 _ v* - ~y^~y^~y^ + y^ + y^ + y^ + ^H ■ (A35) 

\ ^BOd ^BOo ^COti ^COo ^Bj ^B ^C d ^C / 

We can evaluate the properties of the competing phases by determining the minima of the potential Q. In the 
thermodynamic limit, where M is taken to infinity, these minima will give the exact solutions for the system since a 
saddle-point evaluation gives limM_ >00 (l/8A/)lnZ = — min^A^)- 

We mentioned earlier an alternative to the bipartite symmetry This alternative choice consists in taking the matrix 
elements between states pi + Q and P2 + Q equal to those between pi and p2, so that Eq. (|2"2")) becomes 

H„ = ( D E f % ) , (A36) 
where now are complex. Such a choice modifies the four-fermion potentials Yboo, Ycoo,Yb o , and Yc a as 

Ym ° = 8ME 2 E( 1 ' 2 J+ 4 J 3 ')( h - c -)> ( A3? ) 

BOo ij 

Ycoo = ^(li4 j + 2 i 3 i )(h.c), (A38) 

CQo ij 

Yb " = ^^E( 1 « (T2 ^+ 4 J' T3i )-( h - c -)' ( A39 ) 

Bo ^ 

Yc ° = 8M^? (li<T4j ' + 2i °' 3i) ' (h - C - ) - (A40) 
The coupling constants oaf, flsc-s, and asc-d become 

_ 1 ( 1 1 + 1 + 1 1 ( A41 ) 

256M y E sod £ cod E Bd E Cd J ' 

1 / 1 + 2 + 1 + 2 3 6 3 6 \ ^ A42 ^ 

256A/ yS| orf S| 0o ^C0d ^COo ^B d ^B ^Ci ^c / 

1 / _L 2_ + _L 2 3_ + _6 3_ + _6_\ ^ A43 ^ 

256M yS Bod S| 0o S cod S^, 0o S 2 ^ E 2 -^ £c / 

Again, the interaction can be made attractive in the AF channel by favoring spin over density exchanges, while it its 
attractive in the <i-wave channel and repulsive for s-wavc pairing when large momentum transfers are favored. The 
resulting thermodynamic potential has the form of Eq. (|A33[) . although with slightly different expressions for A and 
B. The main results in the text remain valid with this alternative choice. 
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